Proceedings of the 6th European VLBI Network Symposium 

Ros, E. & Zensus, J. A. (eds.) 

June 25th-28th 2002, Bonn, Germany 




EUROPEAN 
NETWORK 



Higher Resolution VLBI Imaging with Fast Frequency Switching 

E. Middelberg,^ A. L. Roy/ R. C. Walker,^ H. Falcke^ and T. P. Krichbaum^ 



O 

o 

(N 



00 
O 

o 

O 



Oh 
I 

O 



^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
^ National Radio Astronomy Observatory, P.O. Box 0, Socorro, NM 87801, USA 

Abstract. Millimetre- VLBI is an important tool in AGN astrophysics, but it is limited by short atmospheric 
coherence times and poor receiver and antenna performance. We demonstrate a new kind of phase referencing for 
the VLBA, enabling us to increase the sensitivity in mm- VLBI by an order of magnitude. If a source is observed 
in short cycles between the target frequency, Ut, and a reference frequency, v,-c{, the v% data can be calibrated 
using scaled-up phase solutions from self-calibration at i/ref . We have demonstrated the phase transfer on 3C 279, 
where we were able to make an 86 GHz image with 90 % coherence compared to self-calibration at Vt- We have 
detected M81, our science target in this project, at 86 GHz using the same technique. We describe scheduling 
strategy and data reduction. The main impacts of fast frequency switching are the ability to image some of the 
nearest, but relatively weak AGN cores with unprecedented high angular resolution and to phase-reference the Vt 
data to the v^^f core position, enabling the detection of possible core shifts in jets due to optical depth effects. 
This ability will yield important constraints on jet properties and might be able to discriminate between the two 
competing emission models of Blandford-Konigl jets and spherical advection-dominated accretion flows (ADAFs) 
in low-luminosity AGNs. 



1. Introduction 

The regions where AGN jets are launched and coUimated 
are difficult to observe with VLBI because the jets are 
launched very close to the black hole and most of the 
bright objects are very distant. Only in the closest AGN 
in Sgr A*, M87, M84, Cen A and M81 can the highest res- 
olution observations resolve several tens of Schwarzschild 
radii, comparable to the scale of 10-1000 Rs where jets are 
predicted to be launched and coUimated (eg Koide et al. 
2000, Appl & Camenzind 1993). Whilst the coUimation 
regions remain unresolved in Sgr A* and M81, there are 
hints for resolving it in M87 (Junor et al. 1999). 

The highest resolution VLBI images are currently 
achieved from observations at 86 GHz. However, short at- 
mospheric coherence times and poor receiver performance 
and antenna efficiencies severely limit the observations to 
sources with S'sgghz 0.4 Jy. We have developed a new 
phase-referencing strategy for the VLBA that is not lim- 
ited by self-calibration at the target frequency vt , but only 
at a lower reference frequency v-ccf- A source is observed 
switching between these two frequencies on a timescale 
that does not exceed half the atmospheric coherence time. 
After self-calibrating the fjcf data, the phase solutions are 
multiplied by the frequency ratio r = vt/v-^ci and added 
to an instrumental, antenna-based phase offset A0. Using 
these phases, the vt data can be imaged. 

In a pilot project, we have observed M81 and 3C 279 
on January 5, 2002, with various frequency pairs and cycle 
times to develop the calibration technique. We describe 
the details and results of this project and what we have 
learnt for future observations. 



2. Observing Strategy and Frequencies 

The Very Long Baseline Array offers the opportunity 
to switch between frequencies in less than ten seconds. 
Frequency switching is done by moving the subreflector, 
and therefore the switching time is limited by its speed 
(the setup of the electronics only takes one or two sec- 
onds). The frequency selection is governed by three con- 
siderations. 1) The feed horns should lie as close as pos- 
sible on the antenna feed circle in the secondary focus to 
minimize the switching times and hence the data loss. 2) 
The phase noise at vt should be as low as possible, and 
one should therefore select v-^ci such that the combination 
of antenna performance (antenna efficiency and receiver 
noise) and source flux density at z/rcf yields little phase 
noise after multiplication by r. 3) The target source needs 
to be strong enough for self-calibration during the envis- 
aged integration time at v-^cf- One should at least expect 
a 5 (7 detection within this integration time. 

We assumed a nominal coherence time at 86 GHz of 
30 s when planning the experiment and therefore selected 
an integration time of 15 s per frequency. We expected a 
data loss of 5 s due to frequency switching, leaving 10 s of 
data at each frequency. Within 10 s, and using 256 Mbps 
recording for higher sensitivity, we calculated the theoreti- 
cal SNRs for three potential reference frequencies 15 GHz, 
22 GHz and 43 GHz (Table |l]). However, the 22 GHz feed 
horn has a large angular separation from the 86 GHz feed 
horn, requiring long switching times, so we did not con- 
sider it further. The phase noise after scaling by r is about 
equal for 15 and 43 GHz, and although the switching time 
from 43 to 86 GHz is shorter, we preferred the higher de- 
tection SNR at 15 GHz. 
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Table 1. Possible reference frequencies, the angular sepa- 
rations of their feed horns from the 86 GHz horn, SEFDs, 
the expected M81 flux density as calculated from an esti- 
mated 150 mJy at 8.4 GHz (Bietenholz et al. 2000) and 
assuming a spectral index of a = 0.2 (Renter & Lesch 
1996), the expected SNR in a 10 s observation, the fre- 
quency ratio to 86 GHz and the expected phase noise after 
scaling by r. 



3. Observations 

Observations were carried out on January 5, 2002, using 
all antennas of the Very Long Baseline Array in excel- 
lent weather conditions at all stations, resulting in atmo- 
spheric coherence times of more than 100 s at 86 GHz. 
We observed 3C 279 and M81 using all possible frequency 
pairs of 15, 43 and 86 GHz, although we did not equally 
distribute the observing time among those pairs. M81 was 
predominantly observed switching between 15 and 86 GHz 
because it is at the limit for self-calibration at 43 GHz. 
The cycle times on 3C 279 varied from 30 s (15 s at each 
frequency) to 90 s, and the cycle times on M81 were 30 s 
throughout. 

4. Data Reduction 

Data reduction was carried out in Bonn, Germany and 
in Socorro, USA to exchange experience with the NRAO 
staff. After loading the data into AIPS, the 3C 279 data 
were inspected to develop a flagging scheme based on the 
frequency switching times. We found that the first 6 to 
7 s of each scan needed flagging, independent of t'lcf . We 
noticed that the subreflectors do not always come into 
position, resulting in pure noise during the affected scans. 
The flagging tables generated by the online systems keep 
track of those events. We have merged the online flagging 
tables with our flagging scheme for optimal performance. 
The data were amplitude calibrated using Tgys and gain 
measurements. Fringe-fitting on 3C 279 was performed at 
each frequency on a short section of good data at all an- 
tennas to compensate for single-band delays and intra-IF 
phase offsets. This allowed averaging across the band dur- 
ing self-calibration in the later stages of data reduction. 

4.1. 3C 279 

We developed the calibration strategy using 10 minutes 
of data switching between 43 and 86 GHz (frequency ra- 
tio r = 2.00) every 30 s on 3C 279. Self-calibration was 
used to determine the phase solutions at 43 GHz. The so- 
lutions were exported to a plain text file and read into a 
Python program to multiply the phase solutions by r and 




Fig. 1. 86 GHz image achieved from a 10 minute obser- 
vation of 3C 279, switching between 43 and 86 GHz every 
30 seconds. Self-calibration was done with a 10 minute so- 
lution interval. The peak flux density is 4.17 Jy; contours 
are drawn from -2 to 4 Jy in steps of 1 Jy. 



reformat the data suitable for reading into AIPS. From 
each two adjacent solutions, we calculated the phase rate 
to allow interpolation in AIPS. Although for this special 
frequency pair, r is exactly 2 by good fortune, we soon 
discovered that the non-integer frequency ratio r = 5.61 
chosen for the bulk of our M81 observations caused un- 
predictable phase jumps after multiplying the phase solu- 
tions by r. This can be understood if one pictures a time 
series of phase solutions at z^iof, where the phase under- 
goes a 360° turn between two adjacent scans. This turn 
of phase will have no effect at i/ref, but in case r is a 
non-integer, the phase at Vt will show a step, because 
<j>rci X r mod 360 ^ 0. The problem can only be solved 
if one keeps track of all phase turns at v-rcf over the whole 
observing run. We therefore implemented data plotting 
and interactive solution editing in our software to remove 
or introduce turns of phase in the phase-time series where 
it seemed obvious. After scaling the solutions by r, the 
calibrated visibility phases were constant with time, but 
an instrumental antenna-based offset remained. We mea- 
sured this offset using self-calibration with a solution in- 
terval equal to the length of the observation, in this case 
10 minutes. After this flnal step, the data were imaged. 

The gain in image quality is illustrated in Fig. 
Fig. 1^ shows the data after self-calibration with a solu- 
tion interval of 10 minutes. This has brought the visibil- 
ity phases to zero on average, but the loss due to atmo- 
spheric fluctuations is significant. Fast frequency switch- 
ing (Fig. H) yields an image that shows 90 % of the flux 
density shown in Fig. ^, which has been self-calibrated 
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Fig. 2. The same data as in Fig. |l| but self-calibration was 
done at 43 GHz and the phase solutions were multiplied 
by r and transferred to the 86 GHz data. The peak flux 
density has increased to 9.07 Jy; contours are drawn from 
-2 to 9 Jy in steps of 1 Jy. 



at 86 GHz. Thus, the coherence loss is only 10 %. Fig. ^ 
shows the phase solutions obtained via self-calibration at 
ft and those obtained via self-calibration at v-cct- 

4.2. M81 

The data reduction was done in the same way as for 
3C 279, and self-calibration with a solution interval of 
15 s at 15 GHz yielded about 97 % good solutions. Due 
to the non-integer frequency ratios, it was not possible to 
integrate longer than 35 min, the longest continuous time 
spent cycling between 15 and 86 GHz. Also, only 8 s out of 
30 s delivered good data because we mostly used a short 
symmetric cycle time, so the sensitivity was not sufficient 
after 35 min (i.e. only 9 min of good data at 86 GHz) to 
determine the phase offset between 15 and 86 GHz. As a 
consequence, we were not able to image M81 and we only 
made marginal detections in data plots. (Fig. ^ and ^. 

However, this result is preliminary, and more work is 
in progress. We plan to use more sophisticated methods 
in determining phase solutions at z^j-cf , e.g. to consider the 
ionosphere's electron content, to estimate the tropospheric 
delay and to use polynomial fits to the phase-time series 
for better interpolation results. Every single aspect will 
only have little effect, but together, they might yield more 
reliable phase solutions at 86 GHz. 
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Fig. 3. The same data again, but phases were determined 
using self-calibration at 86 GHz with a solution interval 
of 30 s, equal to the scan length. The peak flux density is 
now 10.4 Jy; contours are drawn from -2 to 10 Jy in steps 
of 1 Jy 

Fringe-fitted and scaled ptiases on 3C279 
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Fig. 4. Comparison of the 86 GHz phases (solid line) and 
the scaled 43 GHz phases (dotted line) on the LA-NL base- 
line. After phase wraps have been taken out, the scaled 
43 GHz phases follow the 86 GHz phases remarkably well. 

5. Conclusions 

We have developed a new phase-referencing strategy to 
make mm- VLBI observations of sources too weak for self- 
calibration. The strategy is based on rapid changes be- 
tween the target frequency and a reference frequency and 
uses scaled-up reference frequency phase solutions to cali- 
brate the target frequency. We have observed 3C 279 as a 
strong test source and M81 as a weak science target. Using 
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Fig. 5. Sample detection of M81 on a short baseline (Kitt 
Peak to Owens Valley) in 20 minutes of data. The top 
panel shows the phase of the vector-averaged cross-power 
spectrum in degrees, the lower panel the amplitude in mJy. 
Eight channels in each IF have been averaged for this fig- 
ure. 
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Fig. 6. Sample detection of M81 on a long baseline 
(Mauna Kea to Owens Valley) within 20 minutes of data. 
Axis are the same as in Fig. ^ 



this technique, we were able to image 3C 279 at 86 GHz 
with a 10 % coherence loss compared to conventional self- 
calibration, and, for the first time, we were able to make 
a marginal detection of M81 at 86 GHz. 

We have learnt three basic lessons from our project: 1) 
To prevent loss of coherence occuring from phase wraps 
at I'rcf, one should select frequencies such that their ratio 
is an integer. Given the VLBA frequency agility, one can 
almost always tune the frequencies to satisfy this condi- 
tion. 2) To increase efficiency, one should use asymmetric 
switching times. The symmetric times that we used gave 
us only 8 s, or 25 % of data per cycle of 30 s duration. 
For future projects, we will use cycles of 60 s with 16 s at 
z/rof , 7 s for switching, 30 s at vt and 7 s for switching back 
to t'rof. This is only possible in good weather conditions, 
when the atmospheric coherence time is 60 s or longer, but 
the VLBA's dynamic scheduling ensures a high probabil- 
ity. 3) One should use regular scans on a strong calibrator 
to monitor the phase offset between the frequencies of in- 
terest. 

Fast frequency switching turns out to be a powerful 
tool for highest resolution VLBA images of weak sources. 



In some of the nearest, but unfortunately radio- weak 
AGNs, it might reveal the jet coUimation region on scales 
of tens of Schwarzschild radii. However, it can also be 
used at lower frequencies, eg. to phase-reference 15 GHz 
observations using interleaved 5 GHz scans. Atmospheric 
coherence and frequency switching times should be less 
problematic at these frequencies. 

Fast frequency switching is basically a phase refer- 
encing technique, so it will reference the vt data to the 
core position at ^'rcf. After phase transfer, the ft phase 
on a given baseline will contain a time-dependent offset 
which comprises the sum of a constant instrumental term, 
A$instr, and a term due to the source geometry at vt, 
A$gco. If the core at i^t is shifted with respect to the 
t'ref core, A'I'gco will be a sinusoid whose period is 24 h, 
whose amplitude is a measure for the core shift between 
the two frequencies and whose zeros give the direction of 
the shift. A<I>instr and A$goo can be separated by moni- 
toring A<I>instr on a calibrator and fitting a sinusoid to the 
residual phase solutions on the source. 

Using this technique, it should be possible to look for 
core shifts due to optical depth effects in radio jets, be- 
cause the bulk of emission in jets comes from the r = 1 
surface, and this surface is expected to move closer to 
the central mass when higher frequencies are observed. 
In low- luminosity AGNs like Sgr A* and M81, two emis- 
sion models have been suggested, using either scaled- 
down jets like in quasars (Falcke et al. 1993), or spherical 
advection-dominated accretion flows (ADAFs; Rees 1982; 
Melia 1994; Narayan et al. 1995). These models can be dis- 
tinguished because only the jet-like emission should have 
a core shift, and we have proposed to look for it in Sgr A* 
and M81 using a variety of frequencies. 
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